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Abstract 

A new method is presented for calculation of the shell correction with the inclusion of the 
continuum part of the spectrum. The smoothing function used has a finite energy range in contrast 
to the Gaussian shape of the Strutinski method. The new method is specially useful for light nuclei 
where the generalized Strutinski procedure can not be applied. 

PACS numbers: 21.10.Pc,21.10.Ma,21.60.Cs 
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I. INTRODUCTION 



Nuclei being far from the bottom of the stability valley are studied extensively at the 
experimental facilities with radioactive beams. One of the fruit of these type of research is 
the production of the light exotic nuclei. Let us refer to e.g. a recently identified new double 
magic nucleus the O l| at the neutron drip line. The exact location of the particle drip 
lines limits the region for these studies and it is intensively investigated both by experimental 
and theoretical methods. Theoretical prediction of the drip lines is based on mass (binding 
energy) calculations since particle separation energies can be easily deduced. 

There are two important theoretical frameworks for global mass calculations. Microscopic 
HF or HFB calculations with sophisticated effective density dependent interactions are very 
successful in this field. In the best HFB mass formula so far [2J the rms error is 674 keV 
j^]. In earlier HF calculations (4, 5] this number was somewhat larger, namely 805 keV and 
822 keV 3]. In order to achieve this improved fit a new parameterization of the effective 
nucleon-nucleon interaction has been introduced and the pairing interaction was treated 
differently than in the earlier calculations. 

Surprisingly a more simple alternative procedure in the framework of the so called macro- 
scopic microscopic (MM) formalism can compete with the microscopic calculations in the 
calculation of the binding energies. The rms error in the MM calculation is 676 keV. We 
may say that the quality of the microscopic and MM methods are the same. Despite the 
almost identical global fits however the microscopic and MM methods show considerable 
differences when the neutron drip line is approached [3j. 

The key quantity of the MM calculations is the shell correction. The concept of the 
shell correction was suggested long time ago by Strutinski y, |7j] and it is still in use. E.g. 
in a recent global mass calculation [8j the basic ingredient of the shell correction method 
the smoothed single particle density is calculated in a semi-classical way by the Wigner- 
Kirkwood expansion. The other elements of the Strutinski method was not altered. 

Since the invention of the shell correction there were several refinements of the original 
method. Besides the original energy averaging, a smoothing in the particle number space 
was introduced 0, [h]]. Even a combination of the two averaging spaces was considered 11]. 
The particle mean field, the simple harmonic oscillator or Nilsson potential was replaced 
in the calculations by more realistic phenomenological forms in which the spectrum has a 
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continuum beside the discrete single particle levels. The treatment of the single particle level 
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13] but an elegant solution 



density due to the continuum was a long standing problem 

was finally reached HQ. 

A large part of the uncertainty due to the proper choice of the technical parameters 
of the smoothing method has been removed by introduction of the generalized Strutinski 



procedure 15|, |l6) , which made it possible to calculate reliable shell correction values for 



medium and heavy nuclei, where the smoothed level density has a long region with linear 
energy dependence. As it will be discussed in Sec. IV., for lighter nuclei the length of the 
linear region is reduced due to the reduction of the number of the occupied shells and the 
increase of the shell gap. For light nuclei the lower and upper ends of the spectrum distort 
linearity, therefore the method is not appropriate for light nuclei. 

The main goal of this work is to develop a new method which is free from this limitation 
and is applicable for the whole nuclear chart, even in the vicinity of the two drip lines. We 
are solving this problem by introducing a finite range smoothing instead of the infinite range 
Gaussian smoothing used in the Strutinski method. 

The paper is organized as follows. In Sec. II. we recapitulate the formalism of the 
calculation of the shell correction. In Sec. III. we describe the standard Strutinski method 
with the plateau condition. In Sec. IV. we do the same with the generalized Strutinski 
procedure, what we want to replace in this work. In Sec. V. we describe the new method 
with finite range smoothing in details. In Sec. VI we apply the new method for several 
nuclei and calculate shell corrections for neutrons and protons. Finally in Sec. VII. we end 
with the main conclusions of the paper. 



II. CALCULATION OF THE BINDING ENERGY BY USING THE SHELL COR- 
RECTION. 

The binding energy of an atomic nucleus composed of A = N + Z nucleons (N neutrons 
and Z protons) B(N, Z) can be calculated in the microscopic-macroscopic model (MM) as 

B(N, Z) = E macr (N, Z) + 5E(N, Z) , (1) 

where E macr (N, Z) is the binding energy calculated in the macroscopic model (e.g. liquid 
drop or droplet model) and 5E(N, Z) is the shell correction. While E macr (N, Z) is a smooth 
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function of the number of nucleons, the shell correction takes care of the shell fluctuations 
of the binding energy which is missing from the macroscopic model. Shell fluctuations are 
present in any microscopic model. E.g. the shell correction can be calculated from sin- 
gle particle energies of self-consistent Hartree-Fock and relativistic mean filed calculations 

nn n 

171 . |18| . In Ref. [18] shell correction calculated on the single particle energies was used 



to generate a smooth energy from the result of these microscopic calculation and the typi- 
cal phenomenological parameterization of the microscopically calculated macroscopic energy 
terms were analyzed. 

In the present work we use the simplest i.e. the independent particle shell model to 
generate the single particle energies in a phenomenological nuclear potential for the sake of 
simplicity only, since the smoothing procedure could be tested equally well on the result of 
this simple model. In this model we treat neutrons and protons separately. In this case the 
shell correction 

6E(N, Z) = 6E AN T ) = 5E(N) + 8E(Z) (2) 

is the sum of the shell corrections SE T (N T ) calculated for neutrons: r = v with N u = N 
and for protons r = 7r with N n = Z . In what follows we shall discuss the calculation of the 
shell correction 5E T (N T ) for a sort of nucleons only. 

The shell correction can be estimated as the difference of the shell model binding energy 
El and its smoothed counterpart E T calculated also in the shell model. 

SE T = El p - E T . (3) 

Here the shell model binding energy 

E : P = zZ E J ( 4 ) 

3=1 

is a sum of the single particle energies EJ of the lowest energy orbits, from E\ until the 
Fermi-level. In the sum above we can take into account the rvfold degeneracies of the shell 
model orbits and use only the different single particle energies denoted by e\ 



Ei P = j2 n * eT i- ( 5 ) 



The key quantity of the MM model is the smoothed energy E T therefore, we have to give 
a unique definition for calculating it unambiguously. If we have the bound single particle 



energies: ej, the density of the bound nuclear levels is 

gl{E) = Y j n l 5{E-el). (6) 

% 

The particle number as a functions of the energy E of the single nucleon considered is an 
integral of the level density in Eq.flS]), i.e. it is equal to the following step function: 



n T (E) 



I s5(e)de = J> - e[) , (7) 

J — oo 



where Q(x) is a Heaviside function of the form: 

f , if x < 

Q(x) = { (8) 

[ 1 , if x > . 

Since in the smoothing procedure we treat neutrons and protons on the same footing, we 
can drop the r index for a moment. ( We shell include it later again when it is needed to 
avoid ambiguity.) We can calculate the smoothed level density g(E) from the level density 
in Eq.flE]) by folding it with a properly selected smoothing function: f p (x). The smoothing 
function spreads the energy of a discrete level over a certain energy range characterized by 
the smoothing range parameter 7. Therefore, the smoothed level density is 

9(E) = ~ [ +0 ° 9(e) f P ( 6 -— -)de. (9) 



7 J -00 V 7 

The smoothing function in Eq.flH]) is usually a product of a weight function w(x) and a 
polynomial h p (x) of degree p 

f p (x) = w(x) h p (x). (10) 

The later is called as curvature correction polynomial. Since the smoothing function f p (x) = 
f p (—x) is an even function of x, for an even weight function w(x) the polynomial h p (x) 
should also be even and the coefficients of the odd terms in it should be equal to zero. 
Therefore, the curvature correction polynomial has the form: 

h p (x) = °i xi ■ ( U ) 

i=0,2,..., P 

The Cj coefficients of the curvature correction polynomial h p (x) are determined from the so 
called self- consistency condition which requires that the smoothing should reproduce 
the original function if it is a polynomial g n (x) with degree n < p + 1: 

/+00 
9n(x') f p (x - x')dx' . (12) 
-00 



We calculate the smoothed energy by using the smoothed level density in Eq.® : 

E = I e~g(e)de. (13) 

J — oo 

The smoothed Fermi-level A is calculated from the condition that the number of neutrons 
and protons, i.e. the particle number is given: 

N = I g(e)de . (14) 

J — oo 

The smoothed Fermi-level A is different from the Fermi-level A because the level density has 
been modified by the smoothing. 

III. STANDARD STRUTINSKI METHOD WITH PLATEAU CONDITION 

Strutinski used a smoothing function with a Gaussian a weight function 

1 

and it can be shown that the curvature correction polynomials for a weight function of 
Gaussian shape are the associated Laguerre-polynomials 

h p (x) = L l J^) . (16) 

Therefore, in the standard Strutinski method the smoothing function is 



w(x) = — 7=exp(— x ) , (15) 



/ P (x) = ^exp(-x 2 )L^(x 2 ) . (17) 



1 

7T 

For nuclei lying on the bottom of the stability valley the single particle potential can be 
approximated by a simple harmonic oscillator (h.o.) form. For a nucleus with mass number 
A the distance of consecutive shells can be expressed by the well known rule {20 1 

hn = 41 A~ 1/3 [MeV] . (18) 

Shell structure of this simple h.o. model is modified by the presence of the spin-orbit 
interaction and also by the non-spherical shape of deformed nuclei but the quantity in Eq. ( ITBl 
is still serves as a reasonably good measure for the shell structure. An attractive feature of 
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the h.o. potential is that the shell correction 8E{^y,p) as a function of the smoothing range 
7 shows a wide plateau in which the 

d5E ^ =Q (19) 
cry 

plateau condition is fulfilled. More precisely, the fulfillment of the plateau condition is valid 
if at the same time the values belonging to the plateau are practically independent of the p 
value used. It was observed that the plateau condition is fulfilled for h.o. potential. Since 
7 and p are technical parameters of the smoothing procedure and they have no physical 
meaning, it is natural to expect that the definition of the smoothed quantities should not 
depend strongly on these values. Therefore, the shell correction calculated for the h.o. 
potential is well defined. This nice feature of the h.o. potential is related to the fact that 
this potential has only bound states (even at high positive energy values). For potentials 
which are similar to the harmonic oscillator potential e.g. the Nilsson potential we can always 
find regions for 7 where the plateau condition is fulfilled 12|, |21[ . Since these potentials have 
only bound states (infinitely many) and no continuum the ending of the bound states does 
not spoil the picture. 



IV. GENERALIZED STRUTINSKI PROCEDURE FOR SPECTRA WITH CON- 
TINUUM 

However a more realistic single particle potential has a discrete spectrum with finite 
number of bound states < and a continuum of scattering states with E > energy. The 
full level density in this case is a sum of the level densities of the discrete states gd{E) and 
that of the scattering states g c {E) forming the continuum 

g(E) = g d (E)+g c (E). (20) 

Now the smooth level density has to be calculated again with the prescription of Eq.([9]). 
It was realized by Brack and Pauli|2l|] that for this case the plateau condition can not be 
satisfied since the 8E(j,p) curves, what we call plateau curves do not have wide plateaus, 
where Eq. (119ft is fulfilled. They searched for the minima 8E{^f pi p) of the plateau curves for 
each p values and introduced the concept of local plateau condition. At the minima i.e. at 
7 = 7 P Eq.( fT9l is certainly satisfied. An additional requirement of the local plateau condition 



7 



is the approximate p-independence of the 5E( r y p ,p) values, which is satisfied if the variation 
of the 8E(j p ,p) values are small. 

It was shown in Ref. that sometimes even the local plateau condition might not be 
fulfilled and the smoothing procedure of the standard Strutinski method might not able to 
furnish us with well defined smoothed energy. A typical nucleus for which the local plateau 
condition fails if the continuum part of the spectrum is taken into account is the 146 Gd, as 
one can see in Fig. HJ Although one can find minima for each plateau curves, the shell 
correction values at these minima vary too much (even an approximate p-independence is 
not hold). Therefore it is not surprising that the 5E(y p ,p) values deviate considerably from 
the semi-classical value. 

In order to cure this difficulty in the work 15| a modified plateau condition was suggested. 



x 2 



In the modified plateau condition the plateau condition in Eq. ffl9|) is replaced by the require- 
ment that in a certain energy region the smoothed level density should be fitted well by a 
straight line. 

The shell correction 5E(j p ,p) for a given p should be calculated with those 7 P value for 
which the smoothed level density can be fitted best by a linear function: y(E) = aE + b in a 
certain energy range: [ej, e u ]. So we should find the minimum of the function in the variable 
7 for each p value 

flu <2 

ii,p) = 7, p) -y(ii) ■ (21) 

i=l 

Here qi for i — 1, .., n u is a mesh of the energy interval [ej, e u ] used, and 7 P is the value where 
the function \ 2 nas its minimum at a given p- value. To get rid of the shell fluctuations the 
length of the interval has to be larger than the estimated shell gap 

e u - ei = 1.5 hVt . (22) 

Having selected the proper 7 P value for a set of p values between p min = 6 and p ma x = 14, 
the mean value and the variation of the corresponding 5E( , y p ,p) values have to be calculated 
as follow: 

SE = (p -l . +2) £ 5E ^> Ph (23) 

\ymax ymin \ ^ ) _ . 

P — Pmin iPmin i , • • • tPmax 



a 



{Pmax Pmin + 2) 

P — Prntn ;Pm%n \^ i • • • iPm> 



{5E{ lpiP ) - 8Ef . (24) 



S 



Since in Ref. 



15[ this variation was reasonably small for most of the nuclei, the mean in 



Eg. ( 123]) was used to define the shell correction and the variation in Eg. (I24j) was considered 
as an uncertainty of the method. The procedure described above was called as a generalized 
Strutinski procedure. 

In order to illustrate the use of the modified plateau condition we present the smoothed 
level densities for the im Gd nucleus in Fig. [2j The lower and upper ends of the energy 
interval in which the best linear fit of the g{E) is required are shown by filled triangles on 
the i?-axis. Practically no p-dependence of the g{E) curves can be observed in the [e/,e u ] 
interval where g(E) is apparently behaves as a linear function of E. Some p-dependence can 
only observed at around E ~ — 10 MeV being a bit above the A value and at higher energy in 
the E = MeV region which has no influence on the shell correction. The large bump of the 
smoothed level density around E — MeV is the effect of the higher end of the spectrum. 
In the positive part of the spectrum only a few neutron resonance contribute to the level 
density and their effect is smoothed by the smoothing parameters which are the abscissas 
of the filled circles in FigJTJ These 7 P values are between 10 — 15 MeV, therefore the end 
effect is spread well below the threshold. The effect of the lower end is less pronounced but 
can be seen at E < —35 MeV. Here the derivative of g(E) with respect to E changes and at 
E < —45 MeV g(E) goes below zero for a while. The main feature of the g(E) is that the 
linearity required in Eq.(JTJ holds only in a certain distance from the lower and upper ends 
of the spectrum. 

In FigJUthe filled circles on the different p curves show the ( i y p ,SE n ( n f p ,p)) points where 
the 7 P values are those where the function in Eq. (1211) has its minimum. One can see from 
the circles that these shell correction values have much smaller variation (a) than the shell 
correction values at the minima of the curves. Moreover the mean of the 5E n ( / ~f p ,p) values 
denoted by circles is in good agreement with the dotted line showing the semi-classical 
value. In the work [15J it was found that this situation is quite typical and the generalized 
Strutinski procedure gave similar values to the result of the semi-classical averaging based on 
the Wigner-Kirkwood expansion in those cases in which the later could be applied. 

Moreover the generalized Strutinski procedure gave similar results to that of the standard 
one for all cases where the plateau condition is fulfilled. But it gave a well defined value for 
the smoothed energy even in cases like 146 Gd where we can not really speak about plateau. 

It turned out only later, in the work [3| where the generalized Strutinski procedure was 
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used for deformed nuclei, that the function in Eg. (1211) might have more than one minimum 
in 7. It was concluded in that the minimum at the smaller 7 value should be selected. 

An uncertainty of the generalized Strutinski smoothing procedure is that the results are 
slightly depend on the position of the [ej,e u ] energy interval used. For medium and heavy 
nuclei the uncertainty of the generalized Strutinski procedure was always below 250 keV. To 
get this small variation, the energy interval [ei, e u ] was adjusted to the smoothed Fermi- level, 
and the upper end of the energy interval was e u = A — HQq. If the interval was shifted up 
to have e u = A and the length was kept the same as in Eq. (l22|) a variation of the shell 
correction by around 400 keV was observed. This uncertainty was still reasonably small and 
it was comparable to the typical deviation from the semi-classical result. 

The dependence on the position of the interval become stronger for light nuclei. If the 
mass number A is reduced, the distance of the shells estimated in Eq. ffl8|) increases and the 
length of the interval in Eq. (l22|) also increases. We should use larger and larger 7 values 
for smoothing the shell fluctuations. On the other hand the region in which g{E) is linear 
becomes shorter and shorter because the effect of the lower end shifts higher and that of the 
higher end shifts lower. Therefore for small A there is not enough space where the required 
linear region could develop. The linearity of g(E) function is spoiled by the end effects. This 
explains why the generalized Strutinski procedure breaks down for light nuclei. 

Therefore, in this work our goal is to find a new smoothing procedure which is less sensitive 
to the end effects, but it still keeps the advantages of the generalized Strutinski procedure 
i.e. the shell correction is practically independent of the p values (a is small). An additional 
requirement is that E resulted by the new procedure should not be too different from the 
result of the semi-classical procedure (Wigner-Kirkwood method) if the later approach can 
be applied. 

V. NEW SMOOTHING PROCEDURE 

A disadvantage of the smoothing procedures used so far is that the Gaussian weight 
function w(x) used has an infinite range, therefore, the effect of an energy e« is smeared to 
the whole energy axis from —00 to 00. Therefore, the effect of the lower and upper ends of 
the spectrum influences the whole region of the smoothed level density and also the shell 
correction 5E. In this work we try to reduce the end effects in these quantities by using 
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weight functions which have only a finite range. One possible candidate for a weight function 
with finite range is a shape 



ke i-* 2 , if \x\ < 1 
w(x) = { (25) 
, if \x\ > 1. 

The value of the normalization constant k should be chosen from the condition that 

1 = J w(x) dx . (26) 

One advantage of the form in Eq. (|25|) is that all derivative of that function are continuous 
at |x| = 1, so the weight function continues smoothly to the regions where it is equal to zero. 
The effect of the smoothing with this form is localized to the x G [—1, 1] interval. In order to 
use the new smoothing function we have to recalculate the curvature correction polynomials 
h p (x) in Eq.l fTTT) for the new weight function (in Eq. ( J25l) ) . The recalculated polynomials 
h p (x) will be different from the one in Eq.( fl6|) and they should satisfy the self-consistency 
condition in Eq. (ll2jl . with the finite-range weight function. As it was shown in Ref. 19J, the 
coefficients Cj of the curvature correction polynomials in Eq. (flT|) are solutions of the system 
of linear equations: 

p 

^2 °i a i+j = ho <j <p , (27) 

i=0 

where the coefficients a; are the integrals: 



a t = J w(x)x l dx . (28) 

The integration is over the interval where the weight function w(x) is different from zero. 

We present the coefficients q for the p G {0,2,4,6} values in Table [I] for illustration 
purposes. In Fig. |3j we present the shape of the smoothing function f p (x) for a few p 
values and the finite range weight function in Eq. (I25[) w(x) = fo(x). In order to show the 
difference to the standard Gaussian case, we present the similar curves with the Gaussian 
weight function in FigJU For both weight functions for p > the curvature correction 
polynomials h p (x) have p = 2m zeroes: 

h p (x ( f) = 0, j = ±1, ±m, x-j = -xj . (29) 

One can observe the positions of the roots Xj of the Eq. (l29j) in Fig. [3] and FigJU For a 
fixed p value it is convenient to arrange the positive roots of Eq.f l29|) so that they form a 
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monotonous series: 

< x { ? < 4 P) < < a# . (30) 

In the smoothing function f p (x) in Eq. ffTO]) the most important part of the smoothing is 
produced by the central region in h p (x): x G [— xf , x^], determined by the first root x^ . 
One can see in the figures that for p > values x^ p+2 ^ < x± i.e. the value of decreases 
when p increases. 

The finite range smoothing has the advantage that the effect of a certain single particle 
energy e$ vanishes beyond the interval E G [e* — 7, + 7]. Therefore, the smoothed level 
density becomes exactly zero for energies lying below (e% — 7), while the Gaussian oscillates 
around zero. This oscillation character appears at any value of the smoothing parameter. 

If we go to higher E- values, we can smooth the oscillatory character of the g(E) if we use 
large enough 7 values in the smoothing function with Gaussian weight function. This is not 
the case however, if we smooth with finite range weight function, where some undulation 
in g(E) remains even if we use large smoothing range parameters. Therefore, it can not be 
well approximated by a straight line as it was in the generalized Strutinski procedure. 

This seems to be an important difference between the smoothed level densities calculated 
by using Gaussian or finite range smoothings. 

We calculate the smoothed energy in Eq. ( 1T31) by using the finite range smoothing 
functions, for a range of 7 G [7 min ,7 mai ] and p G {p m in,Pmin + 2,...,p max } values. This 
allows us to study the plateau curves. For p = the plateau curve is an monotonously 
increasing function, therefore, neither the plateau condition in Eq. (fT9|) nor the local plateau 
condition can be applied. (There is no 7 value where the derivative is zero.) This result 
show the necessity of using curvature correction polynomials. 

For p > plateau curves have minima (and maxima) where the plateau condition in 
Eq. ffT9l) is fulfilled locally. However the plateau curves might have several minima and we 
have to find the proper one among those minima. A necessary condition of the smoothing 
is that the smoothed level density should not reflect the shell structure of the single particle 
levels. Therefore, in the smoothing procedure we have to start searching for the minimum 
of 5E(j,p) from a (p-dependent) 7 min value with which the shell structure has already 
disappeared. 

The most important characteristics of the single particle spectrum is the largest gap 
between the occupied levels. Therefore, we have to determine the largest distance between 
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the consecutive occupied levels of the N particles (shell gap) 




This G value is a more accurate measure of the shell structure of the single particle energies 
than the MIq in Eq.f jT8|) . In order to estimate a reasonable 7 m j ra value, we have to determine 
the effective width of the smoothing function with a given p. The effective width corresponds 
to the central peak of h p (x) in the interval x G Since the effective range of the 

smoothing function decreases for increasing p, therefore, for larger p value one should use 
larger 7 values for having the same smoothing effect. In order to compensate this effect, it 
is worthwhile to introduce a renormalized smoothing range as follows: 

r p = x { ?\ , (32) 

in which the p dependence of the smoothing is considerably reduced. 

In order to smooth the fluctuations due to the major shells this T p range should be larger 
than the shell gap T p > G. To achieve this we introduce a factor F > 1, and calculate 
a minimal value for the renormalized range T P)m i n = FG. (We observed that the optimal 
value for the factor F is F — 1.5 — 2 for light and F = 2.5 — 3.5 for heavier nuclei.) Having 
fixed this minimum we search for the first minimum of 8E(^j,p) for 

1 — 1p,min TpJ- (33) 

X 1 

This criteria serves as a guide to select the proper minimum of the plateau curve 5E{^ p) p). 
For most nuclei the plateau curves have multiple minima at 7 p i < 7 P)2 <,••>< l P ,i- The 
number of minima / generally increases when p increases. We observed that for p = 2 we 
have at most two minima, i.e. I = 1 or I = 2 and one of them satisfies the following condition: 

T 2 j = x[ 2) 72 , ~ FG . (34) 

For higher p values the proper minimum should be close to this value since we reduced the p 
dependence considerably by using the renormalized smoothing range. Therefore, we have to 
select the fc-th minimum, for which T p ^ = X\^Jp,k ~ ^2,1 ■ If we select the smoothing range 
according to this criteria then the variation of the corresponding 8E(j p ^,p) values will be 
small. 
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VI. DETAILS OF THE NUMERICAL CALCULATIONS 



We used Saxon- Woods (SW) potential with spin-orbit term. For protons it was comple- 
mented by a Coulomb potential of uniformly charged sphere with diffuse edge. (To have 
this form is necessary for being able to calculate semi-classical results for comparison.) The 
parameters of the potentials were that of the so called universal potential given in Ref. [28]. 
The depth of the central potential for neutrons (r = v) i 3 = 1/2 or for protons (r = n) 
t 3 = -1/2 

r N — Zi 

V T (Z, N) = -V[l- 2 ^3 — ^ — J , (35) 

where k = 0.86, V = 49.6 MeV. The depth of the spin-orbit potential 

V— (36) 



4 \2fic, 

with the reduced mass p of the nucleon and A so = 35(36) for neutrons (protons). The 
diffuseness was a = a so = ac = 0.7 fm the same for all potential terms. The radius 
parameters were r = 1.347 fm, r = tq = 1.275 fm for neutrons and protons, respectively, 
while for the spin-orbit term r so = 1.31(1.32) fm for neutrons (protons). These potential 
parameters might not be optimal for the individual nuclei but give a good general N, Z 
dependence all over the nuclear chart at least for our purpose for testing our method. 

The single particle energies e, of the single particle Hamiltonian were calculated by diag- 
onalizing the matrix of the Hamiltonian in h.o. basis having twenty principal h.o. quanta 
and maximal orbital angular momentum nine. (An increase of the size of the basis did not 
change the results.) The same basis was used for diagonalizing the free Hamiltonian (with- 
out nuclear potential terms) to get the positive energies ef^ needed to include the effect 
of the continuum in the Green's function method described in Ref. 16| in detail. From the 
difference of the smoothed level densities of the spectra of the true and the free Hamiltonians 
the effect of the artificial nucleon gas cancels out and we get the same smoothed continuum 
level density as we could get by smoothing the continuum level density derived from the 
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derivative of the scattering phase shifts 

In Fig. |5]we show the plateau curves for the im Gd nucleus with the finite range smoothing 
and the result of the Wigner-Kirkwood calculation as a reference. The range of the p values 
used in the present work was taken to be the same as in Ref. 15] in order to make comparison 
with those results possible. Using the new method with the finite range smoothing we are 
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able to use the local plateau condition i.e. to choose the 7 P values where the 8E(j, p) curves 
have minimum for all the plateau curves shown. The shell correction values at the minima of 
the curves agree very well (within 500 keV) with the horizontal line representing the result 
of the semi-classical calculation. Since the a variation of the 5E(jp,p) values in Eq. (|24p is 
small the shell correction value calculated from the mean in Eq. (l2Uj) is well defined. 

In Figj6]we show an example for the double magic 132 Sn nucleus where the a variation is 
smaller that 200 keV and the deviation from the semi-classical value A is less than 1 MeV. 
This is the largest deviation from the cases listed in Table HT1 One can observe in both Figs. 
E]and [HI that the 7 P values, where the minima of the 5£ , (7 P ,p) appear are increasing with 
increasing p values. This can be compensated to some extent if we use the renormalized 
smoothing range T p defined in Eq. (lH2l . 

The SEl'jp, p) plateau curves are very similar for most nuclei we calculated if we select the 
values of the first 7 P minima of the different p curves beyond 7 Pim i„ in Eq. (l33|) . We identify 
the shell correction with the mean values of the £.E(7 P ,p) in Eq. (l23|) and its o variation 
with the uncertainty of the shell correction. 

In Table [IT] we show the shell corrections for neutrons and for a set of medium and 
heavy nuclei resulted by the new smoothing procedure 5E n (FR), and that of the generalized 
Strutinski procedure 5E n (G). Their a variations are in the third and in the fifth columns. 
In the last two columns we compare their values to that of the semi-classical procedure given 
in Ref.jisj]. The differences from 5E SC are below 1 MeV for the new procedure which is a 
bit better agreement than it is by using the generalized Strutinski procedure. The average 
of the differences are 0.6 MeV and 0.8 MeV for these two procedure, respectively. 

In Table II III we show the similar results for protons, where the average of the differences 
from the semi-classical results are 0.4 MeV and 0.6 MeV for the new procedure and for 
the generalized Strutinski procedure, respectively. So the new procedure can be applied for 
protons as well. 

These differences are not large neither for neutrons nor for protons. The result of the new 
procedure is generally closer to the semi-classical result if we approach the drip lines. See 
e.g. the 78 Ni, 122 Zr, l24 Zr nuclei for neutrons and the 180 Pb nucleus for proton. Therefore, 
we believe that the finite range smoothing allows us to approach the drip line closer than 
we can approach it by using the infinite range Gaussian weight function. 

The basic advantage of the new method is however, that the determination of the proper 
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shell correction value is better defined. The values resulted by the new procedure are free 
from most of the uncertainties of the generalized Strutinski smoothing procedure. E.g. they 
do not depend on the position of the interval where the linearity of the smoothed level 
density is required. 

The most important advantage of the new procedure is that it can be applied for light 
nuclei where, as we have discussed in Sec. IV. the generalized Strutinski procedure can not 
be applied. 

The results of the new method for light nuclei are shown in Table IIVI for neutrons and 
in Table |V] for protons. One can see that the agreement with the semi-classical values are 
as good it was for heavier nuclei. We received specially good agreement for oxygen isotopes 
even at the neutron drip line. 

In FigJT] we show the neutron plateau curves for the new double magic nucleus 24 as 
functions of the renormalized smoothing range parameter T p , for p = 6, 8, 14. The semi- 
classical result is the dotted horizontal line. The minima of each curve are denoted by filled 
circles on the corresponding curves. One can see that the SE n {T p ,p) values denoted by 
circles are between -0.9 and -2.3 MeV and their T p values are quite similar at T p ~ 8MeV. 
The variation of the 5E n (T p ,p) values are a ~ O.bMeV and their mean value coincide with 
the semi-classical value. This is certainly an accident but one can see that the A value is 
small for the other O isotopes too. Observe also that the positions of the minima of the 
different p curves in this figure scatter much less in F (~ 15 %) than the locations of the 
minima in Figj6] where the smoothing range 7 was used (~ 90 %) or in Figj5] where the 
smoothing range 7 was used (~ 70 %). 

Therefore, we believe that the finite range smoothing allows us to approach the drip line 
closer than we can approach it by using the infinite range Gaussian weight function. 

VII. CONCLUSION 

The new method uses a finite range smoothing function which makes it possible to localize 
the effect of a single particle state with energy to a finite energy range: [e, — 7, + 7]. 
This localization makes possible to extend the region of applicability of the method to closer 
to the end regions of the spectrum. This helps in calculating shell corrections for slightly 
bound nuclei lying closer to drip lines and also for lighter nuclei, where the shell gap is large, 
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therefore, larger values of 7 values are needed to smooth the shell structure out. The new 
method works equally well for calculating neutron and proton shell corrections. 

We introduced a renormalized smoothing range in which the p dependence of the smooth- 
ing range was reduced considerably. Using this renormalized range the selection of the proper 
minimum of the plateau curves was easier. 

Therefore, we recommend the use of the new procedure with finite range smoothing first 
of all for light nuclei, where the generalized Strutinski method can not be applied. We also 
recommend its use in regions being close to drip lines where the finite range smoothing seems 
to work somewhat better than the generalized Strutinski method. 
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Nucleus 5E n (FR) a SE n (G) a 5E SC A FR A G 



68 m 


0.16 


0.12 


0.50 


0.07 


0.81 


0.65 


0.31 


78 Ni 


-3.59 


0.07 


-2.78 


0.16 


-4.21 


0.62 


1.43 


90 Zr 


-7.42 


0.06 


-7.35 


0.17 


-6.82 


0.60 


0.53 


l22 Zr 


-5.92 


0.11 


-4.52 


0.15 


-6.33 


0.41 


1.81 


124 Zr 


-4.12 


0.12 


-3.25 


0.13 


-4.35 


0.23 


1.10 


ioo Sn 


-8.16 


0.20 


-6.95 


0.23 


-7.50 


0.66 


0.55 


132 Sn 


-9.85 


0.14 


-8.58 


0.10 


-8.87 


0.98 


0.29 


li6 Gd 


-10.26 


0.07 


-10.33 


0.20 


-9.79 


0.47 


0.54 



TABLE II: Neutron shell corrections 5E n and their variations a calculated using the finite range 
weight function (FR) and the generalized Strutinski procedure G in comparison with the semi- 
classical shell correction 5E SC = E sc — E™ calculated for several nuclei. Their deviations from the 
semi-classical result Apr = \5E SC — SE n (FR)\, Ac = \SE SC — SE n (G)\ are also shown. All energies 
are in MeV units. 



Nucleus SE p (FR) 


a 


SE p (G) 


a 


SE SC 




A G 


90 Zr 


1.59 


0.19 


1.88 


0.20 


1.44 


0.15 


0.44 


ioo Sn 


-7.47 


0.064 


-7.42 


0.14 


-7.01 


0.46 


0.41 


132 Sn 


-7.39 


0.068 


-6.04 


0.12 


-6.64 


0.75 


0.60 


146 Gd 


4.89 


0.10 


5.28 


0.24 


4.52 


0.37 


0.76 


180 Pb 


-8.94 


0.15 


-7.78 


0.04 


-8.62 


0.32 


0.84 


208 Pb 


-7.57 


0.07 


-6.73 


0.03 


-7.29 


0.28 


0.56 



TABLE III: Proton shell corrections SE p and their variations a calculated using the finite range 
weight function (FR) and the generalized Strutinski procedure G in comparison with the semi- 
classical shell correction 5E SC = E sc — E™ calculated for several nuclei. Their deviations from the 
semi-classical result Afr = \5E sc — SE p (FR)\, Ag = \SE SC — SE P (G)\ are also shown. All energies 
are in MeV units. 
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FIG. 1: Neutron shell correction 5E n ('y,p) for the nucleus 146 Gd as a function of the smoothing 
range 7 calculated for p = 6, .., 14 by using the Gaussian weight function for the smoothing functions 
f p . Filled circles on the different curves denote the (j p , 5E n {^ p )) points, where 7 P values belong to 
the minima of the function in Eq.Q21j) and the 8E n (p/ p ,p) values are the results of the generalized 
Strutinski procedure. Dotted horizontal line shows the value of the semi-classical value 5E SC = 
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E [MeV] 

FIG. 2: Energy dependence of the smoothed level densities calculated in the generalized Strutinski 
procedure for p = 6, 10, 14 by using a Gaussian weight function for the smoothing functions f p for 
the nucleus 146 Gd . The lower and upper ends of the interval [ej, e u ] in which the condition of the 
best linear fit is applied are shown by triangles on the E— axis. 
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FIG. 3: Shapes of the finite range smoothing function f p (x) for p = 0, 2, 4, 14. Note that fo(x) 
w{x). 



Nucleus 5E n a 5E SC A 

16 -1.63 0.04 -1.57 0.06 

ls O 2.67 0.04 3.01 0.34 

20 O 3.25 0.24 3.11 0.14 

22 C> 0.12 0.53 0.09 0.03 

24 -1.68 0.49 -1.69 0.01 

20 Ne 3.07 0.56 3.01 0.06 

40 Ca -1.77 0.35 -0.66 0.97 

48 Ca -2.91 0.24 -2.59 0.32 



TABLE IV: Shell correction 5E n , the variation a in Eq. (|24p and the semi-classical shell correction 
5E SC = E sc — Eg calculated for several nuclei. The deviations A = \E SC — E\ are also shown. All 
energies are in MeV units. 
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Nucleus 


SE p 


a 


SE SC 


A 


16q 


-1.65 


0.03 


-1.44 


0.21 


18 G 


-1.65 


0.10 


-1.66 


0.01 


20 O 


-2.09 


0.19 


-1.90 


0.19 


22 G 


-2.30 


0.15 


-2.14 


0.16 


24 G 


-3.10 


0.66 


-2.36 


0.74 


40 Ca 


-1.62 


0.12 


-0.91 


0.71 


48 Ca 


-1.70 


0.19 


-1.44 


0.26 


48 iVi 


-0.80 


0.36 


-1.23 


0.43 


56 iVi 


-3.67 


0.29 


-3.45 


0.22 



TABLE V: Shell correction 5E p , the variation a in Eq. (|24p and the semi-classical shell correction 
5E SC = E sc — Eg P calculated for several nuclei. All energies are in MeV units. 
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FIG. 5: Neutron shell corrections 8E n {p/,p) for the nucleus 146 G<i as a function of the smoothing 
range 7 calculated for p = 6, ..,14 by using the finite-range weight function for the smoothing 
functions f p . Dotted horizontal line shows the value of the semi-classical value 5E SC = E sc — E r s l p . 
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FIG. 6: Neutron shell corrections 8E n {p/,p) for the nucleus 132 Sn as a function of the smoothing 
range 7 calculated for a set of p- values by using the finite-range smoothing function f p . Dotted 
horizontal line shows the value of the semi-classical value 5E SC = E sc — E™. 
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FIG. 7: Neutron shell correction 5E n (T p ) for the nucleus 24 as a function of the renormalized 
smoothing range T p calculated for a set of p- values by using the finite-range smoothing function 
f p . Dotted horizontal line shows the semi-classical value: SE SC = E sc — E" ' . 
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